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STEADYSTATE  ELECTRON  DENSITIES  FOR  IRRADIATED  PLASMAS 


ABSTRACT 

Steadystate  electron  density  profiles  have  been  calculated  for  dif¬ 
fusion-  and  recombination-controlled  plasmas  irradiated  by  an  external 
source  of  electromagnetic  flux.  In  the  former  case  the  internal  electron 
generation  caused  by  progressive  beam  absorption  is  balanced  by  electron 
diffusion  across  the  plasma  boundaries.  In  the  second  model  the  beam- 
stimulated  electron  growth  is  limited  by  recombination. 

OJ) 

The  efficiency  of  electron  diffusion  and  consequent  wall  recombina¬ 
tion  as  a  heat  transfer  process  is  calculated  and  found  never  to  exceed 
307..  Furthermore,  the  nonlinear  diffusion  model  gives  rise  to  a  stability 
criterion  below  which  no  steadystate  electron  configuration  is  possible. 

(U) 


1 


I .  INTRODUCTION 


(1  2) 

Previous  reports  ’  have  dealt  in  detail  with  the  Dynamic  Shield 
Response  (DSR)  of  a  slightly  ionized  plasma  when  subjected  to  electro¬ 
magnetic  radiation.  The  general  features  of  the  interaction  are; 

(1)  a  strong,  initial  thermal  coupling  just  anterior  to  the  resonant 
plasma  frequency  depth, 

(2)  a  local  increase  in  electron  density  coincident  with  the  heat 
deposition  which  in  turn  affects  the  pattern  of  radiation  attenuation, 
and 

(3)  a  consequent  upstream  movement  of  an  ionization  front  resulting 
from  the  electron  growth  which  inhibits  deep  plasma  penetration  of  the 
impressed  radiation. 

(U) 

These  DSR  characteristics  are  displayed  in  Figure  1  for  a  nonuniform 
plasma  with  an  initial  exponential  electron  density  distribution. 

(U) 

Models  previously  treated  considered  only  electron  growth  effects, 
i.e,,  all  relaxation  mechanisms  such  as  diffusion,  recombination,  and 
radiation  were  neglected.  In  this  report  the  influence  of  diffusion  in 
altering  the  electron  density  profile  will  be  considered.  Passing  ref¬ 
erence  will  be  given  to  recombination  and  a  time-dependent  growth  limit¬ 
ing  process  as  well. 

(U) 
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Figure  1.  Time  Response  of  the  Electron  Density  of  an  Irradiated 

Plasma. 


II .  THEORY 


By  analogy  with  heat  flow,  we  may  add  to  the  fundamental  electron 
growth  equation  a  diffusion  term  proportional  to  the  divergence  of  the 
electron  density  gradient.  Thus  in  our  onedimensional  model  we  have 

(1) 


where  n  is  the  electron  density 

f  is  the  ionizing  collision  probability 

^  is  the  molecule  ionization  energy 
2 

aE  is  the  electromagnetic  flux  extinction 
D  is  the  diffusion  coefficient. 

(U) 

In  this  first  order  treatment  we  assume  the  diffusion  coefficient 
is  constant,  leaving  its  exact  magnitude  unspecified. 

(U) 

We  seek  steadystate  solutions  of  Equation  (1)  which,  of  course, 
implies  that  3n/^t  =  0.  This  simplifies  the  hyperbolic  partial  differ 
ential  equation  into  an  elliptic  type  involving  only  one  independent 
variable 


(2) 


(U) 
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In  our  model  electrons  are  generated  at  varying  rates  inside  the 
plasma.  The  solutions  of  Equation  (2)  will  give  us  those  configurations 
for  which  the  outward  diffusion  of  electrons  at  the  boundaries  exactly 
balances  their  internal  growth. 

(U) 

Unfortunately,  the  steadystate  assumption  implies  equilibrium  and 
thus  says  nothing  about  the  time  taken  to  approach  equilibrium.  This 
would  involve  a  complete  solution  of  Equation  (1).  However,  some  ap¬ 
preciation  of  the  magnitude  of  the  growth  constant  may  be  gained  from 

(2) 

earlier  nondiffusion  work  '  where  the  time  constant  T  was  shown  to  be 


_^-l  00  2jt  V 

T^X.  0) 


(3) 


where:  F  is  the  flux  incident  from  large  x 
00  ® 


2  2  2 

n^  m  (€^m/e  ) (o)  +v  )  is  the  critical  electron  density 

e,m  are  the  charge  and  mass  of  the  electron 

q  ■  si  (e  /e^)  is  the  index  of  refraction,  the  squareroot  of 
the  ratio  of  the  dielectric  constant  of  the  medium  and 
free  space 

X,a)  are  the  ^  vacuo  wavelength  and  radian  frequency  of  the 
electromagnetic  radiation 

V  is  the  collision  frequency  of  the  electrons  with  the 
numerically  preponderant  neutral  molecule  background. 


(U) 


5 


We  take  for  our  model  an  infinite  plasma  slab  of  thickness  1.  The 
electron  densities  at  the  boundaries  0  and  H  are  to  be  kept  at  zero. 
This  corresponds  to  infinite  sinks  for  the  electrons  diffusing  to  the 
boundaries.  A  monochromatic,  steady  flux  of  frequency  o)  is  incident 
from  left  to  right  on  the  slab  as  indicated  in  Figure  2. 

(U) 


As  in  the  previous  work  we  assume  the  plasma  absorbs  and  transmits, 
but  does  not  reflect.  This  enables  us  to  use  the  WKB  approximation  which 
transforms  Equation  (2)  to 


(4) 


where  K  =  —  is  the  absorption  coefficient,  a  linear  function  of 

qX  03  n^ 

electron  density  in  this  approximation. 

(U) 


Multiplication  by  the  factor  </n  simplifies  Equation  (4)  to 


i  =  2^  V 

A  “  gn  D  nX  03 
c 


where  A  is  a  constant  having  the  dimensions  of  L  . 

(U) 


(5) 
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Our  task  then  is  the  solution  of  the  integrodif ferential  Equation 
(5)  yielding  the  absorption  coefficient  (and  hence  the  equilibrium 
electron  density)  as  a  function  of  position  for  various  plasma  models. 

(U) 


We  begin  with  the  optical  depth  substitution 


which  eliminates  the  integral  but  increases  the  order  of  Equation  (5). 
Thus 


3 

d  £  _  _  exp(-£ )  d£ 

,  3  ~  A  dx 

dx 

(U) 


This  equation  integrates  at  once  to 


exp(-|)  -  exp(-l^) 
A 


(6) 


(7) 


where  the  constant  |  is  the  value  of  |  for  which  K  =  K  ,  i.e.,  (d/c/dx) 
^m  ^  max 

(d^5/dx^)  =  0. 

(U) 

The  integrating  factor  2d5/dx  enables  us  to  perform  the  second  inte¬ 
gration 


8 


2 


i  [exp(-l)  -  exp(-^J]  ^ 


d  E  d  /  dE  \ 

dx  dx^  ~ 


(U) 


This  leads  ininediately  to 

"  (dx)  ""  A  -  |exp(-|^)]  (8) 

where  we  have  determined  the  constant  of  integration  by  requiring  that 
K  =  (d^/dx)  =  0  at  the  boundary  5=0. 

(U) 


Equation  (8)  expresses  the  absorption  coefficient  as  a  function  of 
the  optical  depth.  This  is  not  the  form  we  want.  We  wish  to  express  k 
as  a  function  of  position.  Before  continuing,  we  note  that  the  imposition 
of  the  second  boundary  condition  /c  =  0  at  |  =  5^,  where  is  the  total 
optical  thickness  of  the  plasma,  serves  to  determine  as  a  function  of 
1^.  Thus,  from  Equation  (8)  at  x  =  i  we  must  have 

1  -  exp(-5^)  -  5^  exp(-5^)  =  0  (9) 


The  solution  of  this  transcendental  equation  which  we  shall  require  pres¬ 
ently  is  plotted  in  Figure  3. 

(U) 

Returning  to  the  main  problem,  both  the  absorption  coefficient  and 
position  can  be  expressed  as  parametric  functions  of  the  optical  depth. 
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From  the  definition  of  the  optical  depth  we  have 


where  £  is  the  plasma  thickness. 

(U) 

By  using  Equation  (8)  and  Figure  3,  we  obtain  the  definite  integral 
relating  the  plasma  thickness  to  the  total  optical  depth  of  the  plasma 

i  j  j  =  [l-exp(-l)  -  &exp(-|J]'^^^  d|  (10) 

This  quadrature  has  been  performed  numerically  and  is  displayed  in 
Figure  4.  Although  the  integrand  diverges  as  |  approaches  0  and 
simple  analysis  shows  that  the  integral  remains  finite.  By  appropriate 
expansions  of  the  integrand  the  following  limiting  values  can  be  found; 

i  j  ^  =  n  -Jz  for  6  ^  «  1  > 
i  ||  =  +  2  £n  2  for  |^  »  1. 

Note  the  interesting  result  that  £  (2/A)  remains  finite  as  |  approaches 

Jo 

zero. 

(U) 
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Figure  4.  Total  Optical  Depth  vs.  Grouped  Plasma  Parameters. 


Finally,  in  Figure  5  we  have  plots  of  the  absorption  coefficient 
(electron  density)  against  position  for  various  fixed  values  of  the 
plasma  parameters.  The  ordinate  is  obtained  from  the  relation 


X  J2  ^  [l-exp(-0  -  i  exp(-|J]^^^ 


while  the  corresponding  abscissa  follows  by  numerical  Integration  from 

r  [l-exp(-l)  -  I  exp(-|  )]  d| 


-i 


^  [l-exp(-l)  -  I  exp(-5^)]  d| 


(U) 


Our  theory  for  the  diffusion  model  is  now  complete  for,  given  a 
plasma  thickness  Z  and  the  value  of  the  lumped  parameters,  A,  we  can 
uniquely  specify  the  equilibrium  electron  density  distribution.  From 
Figure  5  xre  note  that  the  equilibrium  ensity  is  symmetric  about  the 
midpoint  for  small  total  optical  thickness  |  ,  or  equivalently,  for  low 

Jv 

values  of  As,  for  example,  the  incident  flux  F^  increases, 

the  equilibrium  optical  depth  grows  as  well  with  the  profile  becoming 
skewed  toward  the  incoming  beam.  The  steeper  slope  adjacent  to  the 
incident  flux  increases  the  electron  diffusion  out  at  this  boundary 
relative  to  the  diffusion  at  the  opposite  surface. 

(U) 

Thus,  qualitatively,  by  increasing  the  electromagnetic  flux  we 
augment  the  diffusion  at  the  outer  boundary  with  little  effect  in  the 
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Figure  5.  Diffusion  Equilibrium  Electron  Density  Profiles. 


plasma  interior.  Thus  the  nonlinear  plasma  response  tends  to  thwart  ef¬ 
forts  to  alter  inner  transport  behavior. 

(U) 


These  conclusions  are  rendered  more  quantitative  by  con  idering  the 
efficiency  of  electron  diffusion  as  an  energy  transfer  mechanism.  We 
will  assume  that  the  electrons  reaching  the  inner  surface  at  x  =  i  re¬ 
combine,  yielding  up  all  their  ionization  energy  in  the  process.  The 
ratio  of  electrons  reaching  the  interior  boundary  to  all  those  being 
created  in  the  plasma  is,  from  Equation  (7) 


exp(-lnj)  '  exp(-|^) 

Id^/dxlx^O  "  l-cxp(-  -^) 

(U) 


(11) 


On  the  other  hand  the  fraction  of  the  incident  flux  going  into  the 
production  of  electrons  is 


exp(-|)  d^ 


=  f[l-exp(^)]  (12) 

(U) 

The  efficiency  is  the  ratio  of  the  heating  rate  of  inner  wall  re¬ 
combination  to  the  incident  flux  and  is  given  by  the  product  of  Equations 
(11)  and  (12) 
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eff  =  f[exp(-|^)  -  exp(-g^)] 


(13) 


(U) 


The  ionizing  collision  probability  £  is,  of  course,  less  than  unity 

and  for  electron  temperatures  of  several  electron  volts  is  of  the  order 

(2) 

of  0.1  .  Figure  6  is  a  plot  of  the  efficiency  as  a  function  of  equi¬ 

librium  total  optical  thickness.  The  maximum  attainable  efficiency  for 
f  =  1  is  307.  and  occurs  at  =  1.79.  For  large  optical  thicknesses 
the  efficiency  drops  off  inversely  with  leading  to 


eff 


X  _  f  n/  (2A) 


»  1 


(U) 


As  the  total  power  F^  is  Increased  the  total  energy  transferred 
through  wall  recombination  increases  as  the  square -root  a  result 

which  follows  from  the  expression  for  A.  Thus 


fF 

0^  =  yS  =  const  (fF^) 


1/2 


where  is  the  wall  heating  rate. 


(U) 


(14) 


This  conclusion  most  certainly  overestimates  the  rate  of  increase 
of  wall  transfer  with  higher  flux  because  of  the  neglect  of  reflection. 
From  Equation  (7)  we  see  that  the  slope  of  the  density  profile  at  the 
outer  boundary  Increases  linearly  with  the  incident  flux 
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Efficiency 


Figure  6. 
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dx 


for  »  1 


(15) 


x=0 


l-exp(-l^)  1 
A  A 


(U) 


'i'he  discontinuity  in  slope  at  the  boundary  will  tend  to  reflect  the 
incident  beam  with  increasing  effectiveness  for  increasing  flux.  This 
will  temper  the  growth  factor  in  Equation  (14),  so  that  the  squareroot 
wall  transfer  dependence  on  the  flux  should  be  viewed  as  an  upper  limit. 

(U) 


One  unexpected  effect  merits  discussion.  This  is  the  fact  that 
(see  Figure  4)  the  nondimensional  parameter  i  'J(2/A)  approaches  the  finite 
value  n  -JY  as  the  equilibrium  optical  depth  tends  towards  zero. 
Physically  this  implies  that  there  are  no  stable  equilibrium  configura¬ 
tions  when  the  following  lumped  plastna  parameters  are  less  than  unity 


2 

£  fF 

o 

Tt^gn  D 
c 


i)X  0) 


(U) 


(16) 


Let  us  imagine  a  plasma  being  Irradiated  with  an  increasingly  in¬ 
tense  beam.  Since  the  ionizing  collision  probability  is  a  monotonic 
function  of  electron  temperature  and  hence  of  incident  flux,  there  will 
always  be  a  power  range  satisfying  this  inequality.  Wall  transfer  at 
this  flux  intensity  would  therefore  occur  only  in  some  unstable,  un¬ 
steady  fashion.  We  are  forced  to  conclude  that  no  steady  diffusion 


18 


pattern  can  be  established  until  the  power  level  determined  by  the  in¬ 
equality  is  exceeded. 

(U) 

We  turn  now  to  treat  another  electron  relaxation  mechanism,  viz., 
recombination  within  the  plasma.  The  fundamental  equation  in  this  case, 
neglecting  now  diffusion,  becomes 

dn  f  _2  2  /ITS 

^  ^  crE  -  pn  (17) 

We  have  assumed  the  recombination  rate  is  equal  to  the  product  of  a 
generalized  coefficient  ^  and  the  square  of  the  electron  density,  reflect¬ 
ing  the  fact  that  deionization  must  Involve  collisions  between  charged 
particles . 

(U) 

The  steadystate  configurations  are  then  the  solutions  of  the  equa¬ 
tion 


2 

n 


(18) 


(U) 


By  making  use  of  the  proportionality  of  the  absorption  coefficient 
and  electron  density,  we  can  transform  this  equation  to 


K  =  K  exp 
o 


(19) 
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where 


K 


2n  V 


2  fF 

)  _._q 


is  the  value  of  the  absorption  coefficient  at  the  outer  boundary  x  >  0. 


(U) 


Using  elementary  analysis  we  readily  find  as  solutions  of  Equation 
(19)  the  family  of  hyperbolas 


K 


/C  X  +  1 
o 


(u) 


(20) 


Figure  7  is  a  universal  plot  of  these  recombination-controlled 
equilibrium  electron  density  configurations.  Once  again  since  the  beam 
interaction  is  strongest  at  the  outer  boundary  we  find  an  electron 
maximum  there  having  the  finite  value 


n  = 


fF  9 
o  2n  ^ 

Rn^P  qX  00 


(21) 


This  maximum  electron  density  is,  among  other  things,  directly  propor¬ 
tional  to  the  electromagnetic  flux. 

(U) 

As  a  final  model  we  treat  a  time -dependent  electron  profile  growth 

(2) 

which  perhaps  logically  belongs  in  a  previous  report  .  Ignoring  all 
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Figure  7.  Reconblnaclon  Equilibrium  Electron  Density  Profile 


relaxation  effects,  the  fundamental  grov^th  equation  can  be  written,  say 
for  our  infinite  sink  model  irradiated  in  the  positive  x-direction 


^  f  .,2  f  „ 

^  g  (tE  =  g  /c  exp 


(U) 


(22) 


Strictly  speaking  this  equation  is  invalid  since  it  predicts,  for 
example  at  the  outer  boundary  an  unlimited  exponential  growth  of  electron 
density.  The  correct  form  includes  a  factor  which  accounts  for  the  fact 
that  the  probability  of  an  ionizing  collision  decreases  as  the  medium 
becomes  ionized.  Thus  we  can  more  correctly  write  the  growth  equation 
as 

?  O'  n)  ('/o  '' 


where  N  is  the  number  density  of  neutral  molecules  and  positive  ions 
(assumed  singly  ionized) . 

(U) 

A  closed  form  solution  of  Equation  (23)  appears  impossible  to  ob¬ 
tain.  The  solution  at  the  outer  boundary,  however,  follows  readily  by 
elementary  integration 

n  =  - ^ -  (24) 

■0 


[1  +  exp(-t/T)] 

Vo 
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where  is  the  initial  electron  density  at  x  =  0,  Thus  the  plasma  at 
the  boundary  approaches  asymptotically  with  time  the  finite,  completely- 

•  ionized  density  n  =  N. 

(U) 

* 
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III.  CONCLUSIONS 


The  steadystate  electron  distribution  has  been  calculated  for  a 
confined  plasma  in  which  Che  internal  electron  generation  is  balanced  by 
diffusion  across  the  boundaries.  The  profile  is  asymmetric  with  Che 
steeper  slope  at  Che  outer  boundary  where  the  radiation  makes  contact 
with  the  plasma.  An  increase  in  the  flux  density  increases  the  equilib¬ 
rium  total  optical  depth,  steepening  the  slope  at  the  outer  boundary  but 
causing  little  change  at  Che  interior  surface.  The  efficiency  of  wall 
recombination  as  an  energy  transfer  process  can  never  exceed  30%.  At 
best,  ignoring  reflection,  the  wall  heating  rate  cannot  rise  faster  than 
the  squarerooC  of  Che  incident  flux.  If  the  product  of  the  plasma  thick¬ 
ness,  ionizing  collision  probability,  and  incident  flux  falls  below  a 
certain  critical  value,  no  steadystate  electron  configuration  is  possible 
and  Che  profile  becomes  unstable. 

(U) 

For  a  plasma  in  which  the  electrons  are  internally  formed  through 
impact  ionization  and  lost  through  recombination,  the  equilibrium  pro¬ 
file  is  hyperbolic.  The  electron  density  has  its  finite  maximum  at  Che 
outer  plasma  boundary  trailing  monotonically  downward  with  deeper  plasma 
penetration. 

CO 
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